Interaction Quenches of Fermi Gases 
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It is shown that the jump in the momentum distribution of Fermi gases evolves smoothly for 
small and intermediate times once an interaction between the fermions is suddenly switched on. 
The jump does not vanish abruptly. The loci in momentum space where the jumps occur are those 
of the noninteracting Fermi sea. No relaxation of the Fermi surface geometry takes place. 
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Recently, the interest in nonequilibrium quantum 
physics has risen significantly. This is due to exper- 
imental and theoretical progress in treating quantum 
systems with time dependent parameters. No exhaus- 
tive overview of the field is possible but we mention 
the manipulation of atoms in optical lattices fi!| and 
time-dependent transitions between distinct many-body 
states in particular [2] on the experimental side. Sud- 
den changes due to heating by ultrashort laser pulses 
make time-dependent investigations also possible in solid 
state systems [sj. On the theoretical side, there have 
been important advancements in new techniques, e.g., 
time-dependent density- matrix renormalization and 
numerical renormalization [6|[, nonequilibrium dynamical 
mean- field theory (DMFT) [zl,!!,!!] or forward-backward 
continuous unitary transformations (CUTs) [10]. 

The main issues investigated presently are (i) the de- 
scription of nonequilibrium stationary states and (ii) how 
various quenched physical systems approach such states 
or equilibrium states. The existence of conserved quan- 
tities can imply that equilibrium Gibbs states are not 
reached. The systems may approach generalized Gibbs 
states which include the knowledge of the conserved 
quantities, see e.g. Ref. H. In particular, integrable sys- 
tems with their macroscopic number of conserved quanti- 
ties play a special role which may lead to more subtle sta- 
tionary behavior [l^]. The precise way how the station- 
ary states are approached is investigated intcnsivelyb by 
analytical [H, [lj,[l^ and numerical means [16. ITl. fisj. 

In the present article, we do not focus on thermaliza- 
tion itself but on the regime before, i.e., for short and 
intermediate times after the quench. We consider inter- 
action quenches of Fermi gases at zero temperature. We 
study how the ground state |0) (Fermi sea) of a nonin- 
teracting fermionic systems evolves once an interaction 
is suddenly switched on. This question has been inves- 
tigated in the continuum field theory of a Tomonaga- 
Luttinger (TL) model [H, [H in one dimension (ID). In 
leading second order in the interaction U of the Hub- 
bard model in d > 1 dimension, the behavior of the mo- 
mentum distribution (MD) was elucidated by CUTs 



Very recently, nonequilibrium DMFT was used to study 
the issue for d = oo dimensions 0] . 

The latter studies motivate our investigation in partic- 
ular. The numerical data [9] has been consistently ana- 
lyzed in terms of a finite jump An{t) := liuik^kf- n't{k) — 
limfe^fcp+ nt{k) > in the MD 
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which decreases quickly. If the interaction is large enough 
oscillatory behavior in An{t) is found. This was also seen 
in the numerical analysis of spinless ID fermions 

In their study of prethermalization Moeckel and 
Kehrein also discuss the short and intermediate time be- 
havior and state that the fermions at the Fermi momen- 
tum acquire a finite life due to fourth order processes . 
Thus they presume that the jump collapses immediately 
after the quench. The width Afc, over which the jump 
is broadened, is small of the order of {ppU)^ where pp is 
the density-of-states at the Fermi level. 

In view of these contradictory results we aim at eluci- 
dating the behavior of the MD at small and intermediate 
times. We provide evidence that the MD displays a jump 
without broadening. The quintessential reason is that 
long-range correlations of the quenched system remain 
determined by the system before the quench. 

First, we revisit the ID TL model and extend 
Cazalilla's results [3, [l^ to small and intermediate 
times. A finite jump occurs computed for all times in 
all orders of the interaction f25l| . Recall that the scatter- 
ing processes are particularly strong in ID: They destroy 
the conventional Fermi liquid. Thus the survival of a fi- 
nite MD jump in d = 1 can be seen as evidence for the 
persistence of the MD jump also in higher dimensions 
although an alternative view is to take the ID situation 
as being too special to deduce generic behavior in higher 
dimensions. Hence, in order to support our view that 
the ID behavior is generic for the short and intermedi- 
ate times after a quench we secondly discuss the general 
situation in the Heisenberg picture. 

a. Tomonaga-Luttinger Model ID fermionic models 
with linear dispersion and without Umklapp scatte ring 
can be mapped to free bosonic models, see e.g. Refs. |20, 
[2ll . [2^ For simplicity, we first consider spinless fermions 
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with the bosonized Hamiltonian in momentum space 
H = J2 v\<l\blb, + ^ E UiqMibl^bl + b_^b^) (2) 

97^0 q^Q 

where v is the bare Fermi velocity and U{q) is essen- 
tiaUy the Fourier transform of the density-density in- 
teraction (23j . The noninteracting Hamiltonian Hq is 
the one with U{q) = 0. Contributions at g = 
are left out because they do not matter in the sequel. 
The bosonic creation operator is given in terms of the 
fermionic creation (annihilation) operators (cfc) by 

b\ = iy/27r/ {\q\L) J^k A+q^k where L is the length of 
the periodic system. For the MD around the right Fermi 
point at /cp we need the 1-particle correlation 

G^{x, t) = -i(0|7/'R,(a;, i)^R(0, t)|0) (3) 

where '0R(a^i t) = e^*^^ tpB.{x)e~^*^ is the annihilating real- 
space field operator of the right-movers (subscript L is 
for left-movers) at site x (20l. l23j|. Fourier transform of 
G^^{x,t) to momentum space provides the wanted MD. 
The standard bosonization identity [13, HU, IH, reads 

where U^^ is the Klein factor decreasing the number iVp; 
of right-movers by one. The boso nic field ip is given by 
ip^ix) = -a^ E.„9>o W(k|i)&ge^''"""'^'/' where a € 
{R, L} and ctr, = 1 (ctl = — 1 ), Finally, the convergence 
factor a will be sent to 0-1- . 

For evaluating GB.{x,t) we need the time-dependent 
operator ^Fi(a:,i) which is given by ([4]) with ipn{x) being 
replaced by ipa{x,t) = e^*^ ipa{x)e~^*^ . This operator 
can be computed if iJ in ([2]) is diagonalized by the ap- 
propriate Bogoliubov transform bq — Cbq + Sb^-q with 
C = cosh^q and S = sinhflg. The diagonalized Hamil- 
tonian H — Eq + X^ijT^o "^^lll^lbq is characterized by the 
velocity v. For g-independent 9q the calculation is partic- 
ularly transparent. Applying the Bogoliubov transform 
forward and backward [iCt Il4| we eventually obtain 

(pR{x,t) — C^(p-[i{x — vt) + SC(p\^{x — vt) 

-SCipl{x + vt))- S^ipnix + vt). (5) 

Combining this result with the time-dependent version 
of ([4]) and inserting it in ([3]) yields the expectation 
value of the product of exponentials in ipa"^ [x ± vt) and 
(pa^ {±vt) . Such a product can be evaluated by bringing 
the annihilating operators ip to the right and the creat- 
ing operators (p^ to the left with the help of the Baker- 
Campbell-Hausdorfff formula e^e^ = ^a+b+ia,b]/2 
[A, B] is a number only. The required commutators are 
[(pa{x), ip'l^{x')] ~ —\n{2'Ki{(Ta{x' — x) — ia)/L). The final 
result reads 

Gn{x,t) = f'^. . /(^;,r)^^(^+^) (6) 
ZTr(x + la) 

^ \{x + irf-{2vtf\ 

r2 + (2ft)2 r2+a;2 ' ^' 



with 7 := S'2 = C'(C/(0)2). For q independent interaction 
one has r in ([7]) is given by the convergence factor a 
which spoils the proper limit a 0+. Conventionally, 
this is solved by taking into account that the interaction 
is not completely local but has a range r. In equilibrium 
calculations the most convenient assumption is sinh^ Oq = 
7e~''l''l [20L [23|. In the present nonequilibrium context 
it is more convenient to assume sinh^ 0^(1-1- sinh^ ^g) = 
7(l+7)e^''l''l which allows us to discuss Eq. ^ rigorously 
at a = 0-|- for all times t and distances x. In this case, the 
above sketched derivation must be done for each pair of 
modes 5, —q and the commutators finally imply Ij6l7p . For 
comparison, we remind the reader that the equilibrium 
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FIG. 1: (color online) Time evolution of the momentum 
distribution nt(k) around the right Fermi point fcp > 
(AA; :— k — kp) for the spinless TL-model with 7 = 0.1. 

For t = in (HEl), f{v,r) = 1 holds and the MD re- 
sulting from the Fourier transform of e^''^^ / (27r(a:;-|-i0-|-)) 
is the expected step function of the Fermi sea, cf. Fig. [1] 
For < = 00 we have f{v,r) — jir^ -I- x^) which implies 
a power law n(k) - 1/2 oc Afc47(i+7) without a jump. 
Although a stationary correlation has been reached no 
thermalization has taken place because of the simplicity 
of the model: macroscopic number of conserved quanti- 
ties, too little variation in the dispersion [U, [l^. 

We focus on small and intermediate times for which 
we find that G^{x^€) decreases always like Xjx which 
implies a finite jump at the Fermi vector. The time- 
dependent prcfactor of 1/a; is given by the first fraction 
r2/(r2-|-(2wt)2) in ([7]) to the power 27(1-1-7) which agrees 
for < — > 00 with previous analyses [Tj, [l0| . Hence there 
is a completely smooth decrease of the jump 

An(i) = [rV(r' + (2wt)2)]2T'(i+7)^ (g) 

remaining finite at all finite times. The complete MD at 
various times is obtained by numerical Fourier transfor- 
mation; see Fig.[T] In the estimate rjv k, h/ J = ©(Ims) 



3 



J is the hopping element from site to site which leads to 
the stated time scales for atoms in optical lattices [ll, Q ■ 
The MD nt{k) at given k does not evolve monotonically 
in time but displays oscillations in line with previous per- 
turbative |,lQi] and numerical [17] results. 

Inspecting ^ it is obvious why the 1/x proportion- 
ality of G^^{x,t) does not change in the course of time. 
The operator ip^_{x,t) propagates through space at max- 
imum with speed v. Hence there is no way how for a 
given time the long distance behavior is changed. This 
argument relies on the existence of a maximum speed 
■fmax < oo by which information can travel through the 
system. This phenomenon was called light-cone effect by 
Calabrese and Cardy they used the term for the 
propagation of entangled quasiparticles a quench. Note 
for later discussion that such a maximum speed generally 
exists independently of the system's dimension. 

But the prefactor of the 1 /x correlation changes in time 
in spite of the light-cone effect. This represents a major 
change in the correlations. Inspecting our derivation we 
see that this effect stems from local commutators, i.e., 
from commutators between ipa {x ± vt) and tpl^ {x ^ vt) 
or from the corresponding pair at a; = 0. Thus there 
is a multiplicative renormalization of the matrix element 
which links -ip'' {x,t)n\0) to -0t(a;)R|O). Below we show 
that this is also the generic situation in higher dimen- 
sions. Here we point out that for the behavior over short 
distances and short times the TL model is not special. 
All corrections which might be induced by other terms, 
which are present in more general models, will not change 
the qualitative behavior found here because they can be 
treated for short distances and times perturbatively. 

For completeness we wish to extend the results for the 
spinless TL model to its counterpart with spin where 
charge (c) and spin (s) bosonic operators arise from the 
symmetric and antisymmetric, respectively, combination 
of the t and J, bosonic operators [l^, |2l[ . The Hamilto- 
nian is diagonal in charge and spin bosons. It is charac- 
terized by the charge triple Wc^^cTc and the spin triple 
Vs,rs,"fs- In the extended boson identity ^ the sum of 
the charge </?r,c and the spin (pK^s modes occur multiplied 
by I/V2 due to normalization [20]. Pursuing the same 
manipulations as in the spinless case eventually leads to 



Gnix^t) 



27t{x + ia 



(9) 

On the one hand, many qualitative features are the same 
as in the spinless case, in particular the persistence of a 
finite jump in the MD for all times. It is given by 



An{t) = 



7c(l+7o) 



7s(l+7s) 



(10) 

On the other hand, the appearance of two velocities, two 
length scales and two exponents which are possibly dif- 
ferent leads to a richer phenomenology. We refrain from 
showing explicit results because the MDs look very sim- 
ilar to the ones in Fig. [TJ The difference in the decrease 



with the distance x between Eqs. ^ and ([9]) is hardly 
discernible. Hence spin-charge separation, i.e., the dif- 
ference Vc ^ f s is visible only in high precision measure- 
ments. 

b. General Interacting Fermions Because the scat- 
tering induced by interaction is particularly strong in ID, 
the persistence of the MD jump in ID indicates that it 
should persist in higher dimensions as well. Alternatively, 
one may presume the persistence of the jump to be pe- 
cuhar to ID. To support our view that the persistence is 
generic we will analyze the general equations of motion 
below. Note that in all dimensions a maximum speed 
VmnK < 00 exists at which operators can propagate. Thus 
the long range behavior of correlations cannot change 
within a finite times. 

To be specific we consider the 1-particle correlation 
G(r,i) = -i(Ol-0(r,t)-0t(O,t)lO) where the Heisenberg 
time evolution of the operators is induced by the inter- 
acting Hamiltonian H while ]0) is the Fermi sea of the 
noninteracting Hamiltonian Hq. A spin dependence is 
omitted for simplicity. The Heisenberg equation reads 
9fV'^(r, i) — 'i/;^(r, t)]. The commutation with Hq 
(Liouville operator Lq) will propagate the 1-particle op- 
erator only. The commutation with the interaction (Lj) 
generates a particle-hole (PH) pair. The iteration of Li 
increases the number of particle-hole pairs. Hence the 
structure of the solution is 

?/.t(r,t) =p; + pt(pt^t)^ + pt(pt_e-t)2 + ... (11) 

where (_fft) stands for a created particle (hole). It 
is understood that each term in (jlip is normal-ordered 
relative to the Fermi sea ]0). The application Lq will re- 
produce the structure of a term P^P^ H'')'^ , i.e., m stays 
fixed, while Li can increase m by one, leave it unchanged, 
decrease it by one or by two. The first two terms in pip 
are denoted explicitly in d dimensions 



|ri|<i>„,axt 



(12) 



P\P^H^\^ /// /ii(ri,r2,r3,t)- 

J J J\rj\<v„ia^t 

: ^^(ri + r)V3t(r2 + r)^A(r3 + r) : d'^rid'^r2d'^r3. (13) 

The Heisenberg equation generates a hierarchy of cou- 
pled differential equations for the hm{{^j},t). Because 



hm{{rj},0) 



for m > and a term with m PH 



pairs requires at least m applications of Li we know 
hjn — 0{{Ut)™') where J7 is a generic value of the in- 
teraction. The expansion in t can be computed order by 
order. Though we cannot prove the convergence of the 
series t we do not see any means that the radius of conver- 
gence vanishes for Hamiltonians with finite coefficients. 
Certainly, the hierarchy is finite and thus well-behaved if 
calculations up to a finite order in the interaction, e.g., 
t/^, are carried out. 

For G{r,t) the expectation value {0\^J{r,t)li)^O,t)\0) 
must be evaluated which contains terms like 



(0\{HPrP^P^P^H^)i) 



<5™,jG(2"+i)(r,i) (14) 
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because only states with the same number of PH pairs 
can have a finite overlap. The Wick theorem is applica- 
ble because |0) is a Fermi sea so that the many-particle 
correlations can be reduced to products of the initial 
1-particle correlations g{r) G(r,0). Thus we know 
G(2™+i)(r,i) ^ O (g(r)2™+i). The information about 
the jump is encoded in the most slowly decreasing con- 
tribution for |r| — > oo, namely the one for m = 

GW(r,t) = // h*{rut)- 

g{r + n-r2)ho{r2,t)d''nd''r2. (15) 

This double convolution implies n[^\]i) = 
no(k)|/io(k, i)p in momentum space where no (k) g {0,1} 
is the noninteracting MD. Clearly, the nonequilibrium 
1-particle correlations inherit many of their qualitative 
properties from the 1-particle correlations before the 
quench. Interestingly, the jumps in the MD occur at the 
same loci where the noninteracting MD no(k) jumps, 
i.e., at the noninteracting Fermi surface FSq. The 
reduction of the jump is given by 



An{t)\ 



keFSc 



- |/io(k,t)P 



IkeFSo' 



(16) 



The Fourier transform ft.o(k, t) corresponds exactly to hk^ 
in Ref. HO after the forward-backward CUT. These gen- 



eral equations set the stage for the analysis of higher 
dimensional systems and allow us to draw general con- 
clusions. Further analysis has to be numerical and it is 
therefore left to future research. 

c. Conclusions Summarizing, we studied interac- 
tion quenches of noninteracting Fermi gases. The focus 
was the question how the jump in the momentum distri- 
bution (MD) vanishes. In ID the Tomonaga-Luttingcr 
model was investigated quantitatively and in higher di- 
mension the general equations of motions were set up. 
For generic Hamiltonians without diverging coefficients 
we showed that the jump survives for small and inter- 
mediate times, displaying a smooth behavior as function 
of time. If thermalization takes place, we expect that 
the jump decreases exponentially (with or without oscil- 
lations). 

We found that the quenched MD still displays many 
qualitative features of the noninteracting Fermi sea. In 
particular, the loci of the jumps are those of the Fermi 
sea. The Fermi surface of isolated quenched interacting 
fermion models does not evolve at all. The models have 
to be extended in order to incorporate relaxation of the 
Fermi surface geometry. Finally, we point out that the 
approach used here for the general situation can also be 
extended to correlations of two or more particles. 

We gratefully acknowledge many helpful discussions 
with F.B. Anders and S. Kehrein. 
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